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Abstract We establish the long-time asymptotic formula of solutions to 
the (1 + a)-order fractional differential equation qO\ +cx x + a(t)x — 0, t > 0, 
under some simple restrictions on the functional coefficient a(i), where oOl +a 
is one of the fractional differential operators qD^(x'), (qD^x)' = qDI + °x and 
D"(tx' — x). Here, qD" designates the Riemann-Liouville derivative of order 
a G (0, 1). The asymptotic formula reads as [a + 0(1)] ■ x sma u + b ■ xi arge as 
t — > +oo for given a, b G R, where x sma u and xi arge represent the eventually 
small and eventually large solutions that generate the solution space of the 
fractional differential equation qOI + °x = 0, t > 0. 
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1 Introduction 

The present note continues our recent papers [21 El H] devoted to the frac- 
tional calculus variants of several fundamental results from the asymptotic 
integration theory of ordinary differential equations. 

Let us consider the fractional differential equation (FDE) of order 1 + a, 
with a G (0, 1), below 

l O] +a x + a(t)x = 0, t > 0, (1) 

where the functional coefficient a : [0, +oo) — > R is assumed continuous. 
The differential operator qO] +cx is a fractional version of the second order 
operator built by taking into account the decompositions 



x" = (x'Y, tx" = (tx' - x)', t > 



in the ring of smooth functions over (0, +oo). 

To declare the operator ^O l t +a , denote by 1Z£ a ((0, +oo), R) the real lin- 
ear space of all the functions / G C((0, +oo), R) with \im[t 1 ~ a x(t)} G R. 

Recall now the Riemann-Liouville derivative of order a of the function / G 
7e£ Q ((0,+oo),R), namely 



(t - s 



t > 0, 



where T stands for Euler's function Gamma, cf. [51 p. 68]. If the function 
/ is at least absolutely continuous, see [HI p. 35, Lemma 2.2], then the 
derivative exists almost everywhere. Now, we introduce the quantities 

O t =oAo- =j f o D t 



and 



= 0jD q o(f.-- ]<\ R c i. 



The different factorisations [10J of a fractional differential operator might 
lead to some interesting models in mathematical physics. We can mention 
that the fractional differential equations [61 [71 [11] are playing an important 
role in fluid dynamics, traffic model with fractional derivative, measurement 
of viscoelastic material properties, modeling of viscoplasticity, control the- 
ory, economy, nuclear magnetic resonance, mechanics, optics, signal process- 
ing and so on. Basically, the fractional differential equations are used to 
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investigate the dynamics of the complex systems, the models based on these 
derivatives have given superior results as those based on the classical deriva- 
tives, see p. 305], 0E1CE2]. 
Notice that the FDE 

£O l t +a x = 0, t>0, (2) 

has a bidimensional solution space in lZ£ a ((0, +oo), M) generated by the 
smooth functions 1 and t a for i = 1, t a ~ x and t a for % = 2, and t and t a_1 for 
i = 3. 

Regarding the equation (pQ) as perturbation of ([2]), one can ask how close 
a solution a; of ([1]) can get to the solution a ■ x sma ii + b ■ xi arge of (T5]), with 
a, b G R? Some simple restrictions on the functional coefficient a(f) will 
be given next to ensure that an asymptotic formula for the general solution 
of each of the three FDEs exist similarly to the case of classical ordinary 
differential equations. In a loose manner, the formula reads as 

[a + 0(1)] ■ x smaU + b ■ xi arge when t ->• +oo. (3) 

Given the fact that the singular integral operators employed in our proofs 
resemble the integral operators from the two-point boundary value problems 
encountered in the theory of second order differential equations, see pj, we 
think that the Landau symbol 0(1) in our formula cannot be replaced with 
its counterpart o(l) in the majority of circumstances. 



2 The case of ^G] +a 

To establish ([3]), we introduce an integral operator acting within a complete 
metric space and prove that it is a contraction with respect to the space 
metric. The existence of its fixed point will follow then from the Contraction 
Principle and the solution based on the fixed point will obey the asymptotic 
formula. 

We start with a formal derivation of the integral operator. Given x G 
O([0, +oo), R) such that x' G TZC a ((0, +oo), R), we integrate © over [t, +oo) 
to get 

1 f* x'(s) , [ + °° f 
r(l - a) J [t - s) a J t 

where x x = ^m^ ft 0^ds G R. 
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Further, 

i /■* i r x'(ti) 



r(i-a)y (t- s y-°<J ( 8 - u y 



-duds 



— x\ • h / -. rz — / (ax)(u)duds. 

A Fubini-Tonelli argument, see [TTJ p. 29], leads to 

duds 



r(l-a)7 (t-s) l - a J (s-u) 



x'(u) I - ., - . —du 



r(i-a)y J u (t - s y-*( s - u y 

r x '( u ) I T, ds 



B(a, 1 — a 



r(l-a) 



x'(u)du, 



where B is the Beta function, cf. [8j p. 6]. Since B(q,r) = and 
T(l + q) — qT(q), with q, r G (0, 1), we obtain that 

Xi 1 /"* 1 /■ +0 ° 

x(t) = x + — — — • t a + — - / (ax)(r)dTds, 

r(l + a) r(a) y (i - s) 1 J s 

with x(0) = x G E. 

Taking a = Xq, b = f(ffa)' w ^ n a 2 + o 2 > 0, the integral operator reads 



as 



1 I* 1 



T(x)(t) = a + bt a + — ^-^—^J (ax)(r)drds, t > 0. (4) 

Theorem 1 Assume that there exists T > suc/j J!^ 00 s 1+Q |a(s)|d,s < 
+oo and 



max{l,T Q } 
T(l + a) 



- /-T p+oo 

/ |a(s)|ds+ / s a |a(s)|ds 
Jo Jt 



k < 1. 



Taen i/ie FDi? (QJ) /or i = 1 has a solution x G C([0,+oo),M.) wift £ae 
asymptotic formula 

x(t) = a + bt a + 0(t a_1 ) = a + ot a + o(l) waen t -> +oo. (5) 

In particular, 0(1) can be replaced with o(l) m |5|). 
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Proof. Let X be the set of all the functions x G C([0,+oo),R) with 



sup < +oo and d the following metric 



t>T 



d(xi, X2) = max < \\xi — a?2|U°°(o,T), sup 



xi(t)-x 2 (t)\ 



, xi,x 2 g x. 



Obviously, .M = (X, is a complete metric space. 
Notice that 



s 3 \ax\(s)ds < 



s 3 \a(s)\ds + 
C(j)-d(x,0), 



+00 



a(s)\ds 



d(x,0) 



where j G {0, 1}, for every a; G .M . 

Introduce the operator T : M. — > C([0, +00), R) with the formula 
We have the estimates 



< M+T° 



rn+Q! 



and 



\T(x)(t)\ < t c 



\a\+T a \b\ (7(0) 



ri + a 



d(x,0) 



|ax| (s)ds 

te[0,T], 

, t > T, 



which imply that T(x) G .M and 
d(T(s),0) 

< max |l, ^ J (|o| + r°|&|) + max{l, T«} f ^-y • d(x, 0), 

where iGM. 
We also have 

d(T(zi),T(z 2 )) < °^ ,V C(0)-d(gi,aa), x^^gM, 
1 (1 + a) 

which means that T : M. — > M. is a contraction of coefficient 

Let Xq G .M be its fixed point. Following verbatim the computations 
from [31 eqs. (10), (16)], we have the estimates 



(t- s y- a 



\axo\ (r)drds 
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o 

t a 
< — 
a 



*, U ^-{t-TY j ^ /" + °°, i j 

\ax \(T) dr H / \ax \ds 

a a J t 

t 1 f + °° 
\ax \(r) ■ -dr H — / s|ax |(s)(is 



< ^^-d(x , 0) • t a - 1 = 0(t a - 1 ) when t -> +oo. 

a 

Finally, 

x (t) = T(x )(t) = a + bt a + 0{t a ~ l ) when t -»■ +oo. 
The proof is complete. □ 

A particular case of fl5]) has been undertaken in [2] , namely the case when 
a = 1, b = 0. We asked there if, similarly to the circumstances of ordinary 
differential equations [U Section 7] , the solution xq from Theorem [1] would 
have the (powerful) asymptotic behavior 

x (t) = l + o(l) ast^+oo, x' G (L 1 fl L°°)((0, +oo),R). 

We also noticed that, most probably, to get such a result one must look for 
a sign-changing functional coefficient a(t), see [21 Section 3]. 

In the remaining of the present section we shall discuss the issue of 
u x' G L ln and conclude that this can happen (eventually) in very restricted 
conditions. 

Lemma 1 Assume that a G (C fl L°°)([0, +oo),M) verifies the hypotheses 
from fj|/; it has a unique zero t > 0, J °° a(s)ds = 0, J +o ° s\a(s)\ds < +oo 
and B G (L 1 fl L°°)([0, +oo), R), w/jere = t Q ||a|| L oo (ii+oo) /or all t > 0. 
Then, introducing the quantity C{t) = J * jjz^r^ds, t > 0, we have 

\C(t)\dt + sup \C(t)\ < +oo. (6) 
//S* G ^([Oj+ooJjR), w/iere 5*(t) = sup B{s) for all t > 0, then 

s>t 



C*(t)dt < +oo, C*(t) = sup \C(s)l t > 0. (7) 

s>t 

If J + °° s 1+a \a(s)\ds + / +o ° \\B\\ L 2( t +00 )dt < +oo then we also have 

i 

r+oo r+oo / r+oo \ 5 

/ ||C|U 2(U)+oo) du= / / \C(t)\ 2 dt) du< +00. (8) 
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Proof. As in [2], for t > 0, the following estimates are valid 



\C(2t)\<^ + 



a 



as 



o (2* -a) 1 " 



-ds 



(9) 



and 



as 



(2t- s y- a 

t 

a-l 



a(s)ds — (1 — a) 

+ 0O 



a(s)cfe + (l_ a ) y (2f _ a)2 _ a y a(r)drds. (11) 

Since B G L 1 fl L°°, it is obvious that B G L 2 , so we shall focus on the 
second member from the right part of By means of ffTUl) . we get 



a-l 



a(s)ds 



1 



(2t - s) 



a(r)drds 



< t a \\a\\ Lo ° + 



f 2 - 



[s ■ \\a\\L°°)ds = -||a||z/*> ■ t c 



which leads to C G (L 1 n L°°)([0, T ], R), where T = max{l,t }- Further, 
via (fTU), 



a(s) 



(2^-s) 1 - 



-(is 



< t 



< 



a-l 



\a(s)\ds + 1 



a-2 



s|a(s)|ds 



t 2 - 



s\a(s)\ds, t > T , 



and so C G (L 1 nL°°)([T , +oo),R). The estimate flS]) has been obtained. As 
a byproduct, C G L 2 ([0, +oo), R). 

To prove ([7]), introduce D*(i) = sup -D(s) for all t > 0. We rely on the 



estimates 



and 



D*(t)< -WaW^-e, te[0,T o ] 



/ o*W* £2/ £•/ 

-'To </Tn S JTo 



r\a(r)\dT 



To 

^rpa—l p-\-oo 



{I -a) 



T 

r\a(r)\dT, 



To 



Asymptotic integration of (1 + a)-order FDE's 



8 



since the mapping t t— ?■ t a 2 s\a(s)\ds is monotone non-increasing in 
[T ,+oo). 

For the third part, notice that 



and 



9 p+co 

D(t) <-J s 1+a \a(s)\ds, t>T , 



i 

\C(2t)\ 2 dt 

2u 

1 / /' X dt^ * 



< - • \\B\\ L 2 {2Uj+oo) +U -j • 2 s 1+a \a(s)\ds 

< a _1 ||5|| L2(2Ui+oo) + — y s 1+a |a(s)|ds, u > T . 

We have obtained that f+°° (/+°° \C(2t)\ 2 dty du < +oo. 
Finally, 



)2 1 /"To / r+oo \ 2 

du = — / |C(w)| 2 cfo dn 

y/2 Jo \Jau J 

v/2' 

The proof is complete. □ 



< ^^IICIU 2 (o,+oo)- 



Lemma 2 Assume that the function C from Lemma [I] satisfies the restric- 
tions (0) and |2P and either 

||C7|| L » + 2||C7*|Ui = fc 1 <l 



or 



2||C*|| L i < 1, max{||C|U« + \\C\\ L 2, \\C\\ L i + = k 2 < 1, 

where E(t) = ||C|| £ 2( 4j+00 ) /or aZZ t > 0. Then there exists a function y E 
(C fi L 1 n L°°)([0, +oo'),R) swc/i toot 

(/•+oo \ P + CO 

I- J y(s)dsj-J^ (Cy)(s)ds, t > 0. (12) 
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Proof. Set the number 7 > 1 such that 

r+00 

I + 27 / C*{s)ds < 7. 
Jo 

Introduce the set Y of all the functions y G C([0, +00), R) such that 
\y(t)\ < j ■ C*(t), t > 0, and the metric d with the formula 

d(yi,y 2 ) = max{||yi - y 2 |U<»(o,+Do)» Hz/i - ^lU^o.+oo)}, 2/i, 2/2 e V. 

Using the Dominated Convergence Theorem, we deduce that the metric space 
J\f = (Y, d) is complete. 

Consider the integral operator T : A/" — > C([0, +00), R) given by the 
right-hand member of (1121) . The following estimates 



\T(y)(t)\ < \C(t)\(l + \\y\\ L ,) + C*(t) J \y(s)\ds 

< C*(t)(l + 2\\y\\ L i) <C*(t)- fl + 2 7 £ C*(s)ds S j 

< 1-C*(t), t>0, 

show that T : M — > M is well-defined. 
Now, we have 

/+00 
\C{s)\ds-\\ yi -y 2 \\ L ~ 

< {\\c\\ L ~ + \\c\\v)-d{ yi ,y 2 ), 

by noticing that C*(0) = ||C||z,o°(o,+oo)) an d also 

/+00 
(|c(0i-||yi-y2lU0^ 
+00 r+00 

C*(s) I \ yi {r)-y 2 (r)\drds 

r +c 

< 2 

which lead to 



r+00 

/ C*{s)ds-d{y u y 2 ), t>0, 
Jo 



d{T{ yi ),T{y 2 )) < max{||C|U=c + ||C7|| il ,2||C*|Ui}.d(i/i,y 2 ) 
< hd(yi,y 2 ), 



where y\ , y 2 G N '. 
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Notice that we haven't employed (jSj). To do so, let us use different esti- 
mates, namely 



|T(i/i)(t)-T(ito)(f)| < \C(t)\-\\ yi -y 2 \ 



+ 



\C{s)\ 2 ds 



\yi{s) - V2(s)\ 2 ds 



and 



p+oo r+oo 

/ \Vi(s) - y 2 (s)\ 2 ds < sup \yi(r) - y 2 {r)\ ■ \ \y x (s) - y 2 (s)\ds 

Jt T>0 Jt 

< [d(y u y 2 )} 2 , t>0. 

They imply 



\T( yi )(t) - T(y 2 )(t)\ < 



\C(t)\ 



\C{s)\ 2 ds 



d{yi,y 2 ) 



< {\\C\\ L ~ + \\C\\ L >)d(y x ,y 2 ) 



and 



\T(yi)(s) - T(y 2 )(s)\ds < 



\C\\ L i + 



\C{s)\ 2 ds dt 



x d{y 1 ,y 2 ) 



thus leading to 



d(T( yi ),T(y 2 )) < max{||C|Uoo + ||C|| L2! ||C'|Ui + || J B|| Ll }.% 1 ,y 2 ) 
< k 2 d(yi,y 2 ), 

where y x , y 2 G M . 

The operator T '■ M — Y N being a contraction, its fixed point 2/0 is the 
solution of (fT2j) we are looking for. The proof is complete. □ 

Proposition 1 Let y G C([0, +oo), R) be the solution of [W\l from Lemma 
\E If 2/(0) = then the function x G C^QO, +oo),R) with the formula x(t) = 
1 — y(s)ds for all t > is a solution of the FDE (QJ) for i = 1 which 
satisfies the restrictions 



x(*) = l + o(l) as t^+oo, x' G (L fl L°°)([0, +oo), 
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Proof. Following [2], the function x verifies the identity 
If* a(s)x(s) , 1 (* a(s) 

y(t) = -— - / - V \\ J ds = — tzz — r / - ^ ds 

/"+oo 

= -C(*) + / y(r)C(r)dT + C(t) / y(r)dr 
Jo Jt 

= -C(t)(l- £°°y(s)ds}+ J^Cy)(s)ds, t>0. (13) 

We have rescaled C as C{t) = J * , t ^}- a ds, t > 0. 

Let t = in ([12]). Then, = y(0) = - j+°°{Cy){s)ds. This means that 
we can recast the integral expression from f[T3"j) as 

1- y 2/( S )rf S J - J (Cy)(s)ds, t>0, 



which is exactly f|T2|) . 

The proof is complete. □ 

To give some insight to the (still unsettled) issue of u x' G L 1 ", notice that 
the condition y(0) = from Proposition [T] reads as 

r+OQ H air) 

x'(s) / t r- — drds = 0, 



/o Jo 



which is really difficult to handle. A further intricacy is provided by the fact 
that, given a E C([0, +oo),R), the quantity F(t) = f* T^jpU ds, t > 0, does 
not belong to L 1 ([0, +oo), ]R). This follows from 



rt rt p2s \ n ( r \\ /"* Jo r 2T 

L F{2s)ds a M ^h iTds> -Lj^r^-S, HtWt 

— > +oo when t — > +oo, 
where T > is chosen large enough for a to be non-trivial in 2T] . 

3 The case of $O l t +a 

Introduce the relations 

r +oc v(t) 

y(t)=tx'{t)-x{t), x(t) = bt-t ^V^r, t > 0, (14) 



r 2 
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with and y G 7^.£ a ((0, +oo),R), see [I]. 
As before, 



1 



■ / -. r— = Xi + / (ax)(,s)ds, t > 0, 



r(l-a)7 (*-s)< 
where a* = Jim^ /„' G R, and 



y(s)ds 

Xit a 



f + oo 



+ 



r(i + a) r(a) y (* - s 



1-a 



(ax)(r)dTds 



x x t a 1 
+ 



1 



r(i + «) r(a)J (t- s y- a 

r+oo 
X\ + l 



r(l + a) _ 

see HH p. 32, eq. (2.13)]. 
By differentiation, we get 



(ax)(r)dr 



ax)(r)dTds 
' ''" (ax)(v> 




r («) Jo Jo ( s - u) 1 - 



duds. 



y(t) 



j.a-1 



r(a) 



Xi+ (ax)(r)dr 



1 



r* (ax) (5) 

rR7 (t-s) 



1-a 



ds, 



where t > 0. 

Taking a = — ^z^Ffoj an d recalling f JT^|) . our integral operator reads as 



(2-a)r(«) 

T(y)(t) 



t 



a-1 



a + 



r(«) Jo 



j.a—1 r+°° 



T(a) Jo 



ra r 



tar 



sa(s)ds 



sa(s) 



T(a)J (t-s) 1 " 



dudr 



+00 



, , -dudr, t > 0. 

I»y (t-r) 1 ""], « 2 



Theorem 2 Assume that J Q +co t\a(t)\dt + su.pt 1 a J* ffi rfs < +00 and 
1 / /•+" |o(s)|. \ , 



where x — supt 1 a Jq rt-sli-^l 1 -" ^ en ^ e (GP /or z = 3 /ias a 

t>o 1 sj s 

solution x G C 1 ((0, +00), R) with the asymptotic formula 

x(t) = [a + 0(l)}t a - 1 + bt = bt + 0(t a - 1 ) when t -> +00. (15) 
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Proof. Let us start by giving a simple example of x- If the functional 
coefficient a E (C fl i^XfO, +00), M) verifies the restriction 

Ht)\<^, t>o, 



then 



r , m i ds = ( r + n . m i ds 



< 



2f 

dv 



<t'- a 1 } a[s l ds+t 1 -* 1 — -- — ^ 



/- 1 du 



(1 - v) 1 -" 



1 . 2 1_Q 
< - a loo ( 0) i) + a L i(i +00 ) + A < +00, t > 0. 



Notice also that f* ds < t 1 01 f* jj^0^ds < \ and 

Jo (t- s y-<* J (t - s y-« s i-<* 



which leads to the "x" of the mapping t >->■ t 2 ~ a a{t) in [0, +00). 

Introduce now the set Z of all the functions y G C((0, +00), R) such that 
sup t 1_a I < +00 and the metric 



il-a 



d(yi,y 2 ) = supt Q |?/i(t) - y 2 (*)|, 2/1,2/2 e^- 
00 



Observe also that 

"+00 



^r w^M ,, £ 1 . sup( i- lBl(() _, 2(() | 

t>o Jt u z 2 - a t >o 

< d(y 1: y 2 ). (16) 

The metric space V = (Z, d) is complete. Given y e V, we have the 
estimates 

t l - a \T{y){t)\ 

<| a | + J^I [ + °° s \ a ( s )\ ds + ^ -supt 1 - f s ^ s }\ ds 

~ 1 1 r(a)7 1 1 ;| r(«) 4 >s Jo (t- s y-° 
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/ ^^ds ■ sup s 2 ~ a / ^-^du 



w ^ • sup ti- /* - 'y 1 dr ■ sup r 2 - f °° M^! d „, * > 0, 
r(a) t>o 7o (t-r) 1 -^ 1 -" r >o 7 r 



which imply that T(V) Q V . 
Further, we have 



t 1 - a \T(y 1 )(t)-r(y 2 )(t)\ 
< 



' ^la(s)l da+ 1 /■ Kr)| 



r(a) 7o s 1 "" r(a) t >o 7o - r) 1 - a r 1 
X%i,2/ 2 ) 

1 / /" + °° la(s 



r(a) 



^ds + x) % l5 y 2 ), * > 0, 



by means of (TT6I) . where yi, y? £V. 

The operator T : V V being a contraction of coefficient &3, it has a 
fixed point y - Thus, since yo(t) = for large values of t, we conclude 

the validity of the asymptotic expansion f JT5|) for the solution x given by (JUJ). 
Notice also that 

b r +oc 

lim* 1-0 W*) = limt a T(2/ )W = « + ^7-^ / sa(s)ds. 
t\o *\o r(a) J 

The proof is complete. □ 



4 The case of $O l t +a 

The asymptotic formula ([3]) has been already discussed in [3J, however, it is 
worthy to be recalled for reasons of completeness. 

Theorem 3 (TJ1 Theorem 1]) Assume that there exists T > such that 



max{l, T} 
r(l + a) 



^ 1*+^ + / +OC Sa ' a ( S )l rfS ) = ^ < 1 



and J T °° s 1+a \a(s)\ds < +oo. Then, given a, ft 6 1, wi/i a 2 + 5 2 > ; i/ie 
FDi? (fjp /or z = 2 has a solution x G C((0,+oo),R) with the asymptotic 
formula 

x(t) = [a + Oil)}^- 1 + fa Q = bt a + O^ 1 ) w/ien t -> +oo. 
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The formula of the integral operator reads in this case as 

1 ft 1 p+oo 

T(x)(t) = at*- 1 + bt a + — J {t _ s)1 _ a J (ax)(r)drds, t > 0, 

and its fixed point x satisfies also the conditions 

limt 1 - Q x (t) = a, lim { D?x ){t) = T(l + a)b. 

t\0 t— s-+oo 
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